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ABSTRACT 


The time-temperature response of porcine tissue to laser radiation exposure is investigated as a function of wavelength. We 
experimentally measure the thermal response of tissue to laser radiation ranging in wavelength from 1100 nm to 1550 nm. The 
experimental data were compared to simulations performed using thermal modeling software. Based on these simulations, 
and the corresponding experimental data, damage thresholds as a function of wavelength were estimated. This data can be 
used to help optimize the design of optical imaging systems, particularly those being used for biomedical imaging. 


Introduction 

Advancement in laser technology and microscopic techniques, including multi-photon fluorescent microscopy and coherent 
anti-Stokes Raman microscopy,'"^ have provided means to conduct non-invasive, high-resolution, deep-tissue imaging ex¬ 
periments.^ Many of these techniques take advantage of the use of near-infrared (NIR) excitation radiation, which has been 
shown to allow for imaging at millimeter depths in tissue through highly scattering media.Three NIR optical windows 
have been identihed and allow for an increased penetration depth for optical imaging.® The NIR region with wavelengths 
from 650 to 950 nm is most commonly used in optical imaging for biological applications, partially due to the abundance of 
laser sources and detectors that function well in this range.Since Rayleigh scattering, which varies as 1/A^, and Mie 
scattering, which varies as 1/A" with n > 1 for longer wavelengths, decrease with increasing wavelength, there is reduced 
scattering and minimal absorption for NIR wavelengths when compared to visible wavelengths allowing for imaging deeper 
into tissue. The second NIR window has been defined to cover the wavelength range of 1100 to 1350 nm, while the third 
window covers 1600 to 1870 nm. These ranges fall on either side of a strong water absorption band centered near 1450 nm. 

Limited studies have been conducted using the second or third optical windows.'^ With the increasing presence of high- 
powered, short-pulse optical parametric oscillators (OPOs) and optical parametric amplifiers (OPAs) in research laboratories, 
coupled with advances in the spectral response of NIR charge coupled device (CCD) image sensors, additional freedom is 
afforded to researchers when selecting the excitation wavelength for deep-tissue imaging. These windows are of particular 
interest for optical techniques such as coherent anti-Stokes Raman spectroscopy, where exogenous contrast methods are not 
needed for imaging.^ 

The depth at which images can be obtained is partially limited by the amount of excitation radiation that can be delivered to 
the tissue without causing damage. Photothermal damage, as a result of absorption of radiation resulting in localized heating, 
and photomechanical damage, caused by mechanical forces due to the rapid absorption of short pulses of light, are two of 
the primary damage mechanisms for NIR radiation. Damage thresholds at select wavelengths in the second and third NIR 
windows have been published,'^’but further work is required to fully assess damage thresholds over this range. 

In this report, we present the laser-induced tissue heating for porcine tissue in the wavelength range of 1100 to 1550 nm. 
Previously, laser-induced heating from 700 to 1064 nm was studied.'^ The laser-induced tissue heating has been shown to 
be linear with the incident power at a given wavelength;'^ thus, estimates for the amount of power that can be used without 
causing thermal damage can be extrapolated from the experimentally measured thermal response of tissue to a given incident 
energy. While consideration must also be given to photomechanical damage, particularly while using sources with high peak 
power, estimates of the thermal damage threshold are of great value to researchers attempting to push the limits of deep-tissue 
imaging. 



Figure 1 . (A) Rendering of experimental setup used to measure thermal response (some components not shown). OPO, 
optical parametric oscillator; D, dichroic beamsplitter; Camera, FLIR thermal camera; Sample, mounted tissue sample. The 
red beam corresponds to the radiation used for heating (1100 to 1550 nm) and the white beam cotTesponds to the measured 
thermal signal (3 to 5 /im). (B) Sample thermal image taken from a 1425-nm laser exposure shortly after the laser pulse 
atTived. The black circle cotTesponds to the region of interest (ROI) used for calculating a temperature rise. 


Methods 

Experimental setup 

To achieve the desired tunable wavelength range for these experiments, a Coherent MIRA OPO, pumped by a Coherent 
Chameleon titanium-sapphire (Ti:sapph) laser was used as the excitation source. This source produced 120-fs pulses with a 
repetition rate of 80 MHz and a maximum power of 1.2 W over the tunable range of 1100 nm to 1550 nm. The laser light was 
coupled into a IQQ-jim core multimode optical fiber. The distal tip of the fiber was then imaged onto the tissue surface with a 
magnification of 2.5X, producing a 500-/im excitation spot. In order to control the laser power measured at the sample plane, 
the laser was introduced into a A/2 Fresnel rhomb (Thorlabs FR600HM) followed by a polarized beam splitter. The Fresnel 
rhomb allows for A/2 retardance over the broad wavelength range of 600 to 1550 nm, and rotation of this optic allows for 
adjusting the transmitted laser power such that 80±2 mW of average power (42±1 W/cm^) was measured at the sample plane 
for all wavelengths. Fig. 1(A) shows a rendering of the optical system used. 

The tissue sample was imaged to a thermal camera (FLIR SC6000) and thermal prohles were recorded at 200 Hz using 
a 256 by 256 pixel sub-array. The held of view for the thermal images was measured to be 2.4 mm by 2.4 mm with a pixel 
size of 9.3 /J,m. Thermal prohles were recorded in triplicate with 25-nm steps over the 1100 to 1550 nm wavelength range. 
The tissue sample was translated laterally following each exposure to minimize error resulting from residual heating or loss 
of water caused by the previous exposure. 

Tissue samples were acquired using a tissue sharing protocol in place at the Air Force Research Laboratory. Excised 
porcine tissue taken from the Hank of a Yorkshire pig was used as the tissue samples. The tissue samples measured approx¬ 
imately 6 inches square and 1 inch thick, and included epidermis, dermis, muscle and fat layers. A total of 5 samples were 
used in these experiments. Prior to imaging, all samples were placed in an incubator which was maintained at 35 °C with a 
relative humidity of 80%. Each sample was mounted into a custom sample holder. The mount was then placed into the clear 
acrylic box shown in Eig. 1(A) which included a heating element as well as humidity control. The box was maintained at 
35 °C during imaging with the relative humidity inside ranging from 50% to 60%. 

Eor each wavelength examined, the tissue was exposed to the laser radiation for two seconds and the thermal response 
was recorded for a total of 10 seconds. The thermal camera began recording one second prior to laser exposure to provide a 
baseline tissue temperature. Before laser exposure, the tissue sample temperature measured between 29°C to 33°C. 

Thermal simulations 

These experimental results were then compared to simulated thermal prohles using a heat transfer model.'® The BTEC thermal 
model, whose name is derived from its four primary authors, provides a 1-D or 2-D cylindrical coordinate system simulation of 
optical radiation and radio frequency thermal interaction with biological tissue. The code supports the illumination of samples 
by various sources, the temperature response from the linear absorption of optical radiation, and the analysis of subsequent 
damage. Einite difference numerical methods are employed in the solutions of heat transfer. The Crank-Nicholson method 
was employed in the solution for heat transfer. The BTEC simulations can be conhgured with a source term defined by a 
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single or by multiple emitters, such as laser, RF, or broadband sources. The tissue used in these simulations is represented as a 
one-dimensional stack of homogeneous layers along the z-axis.^**’^' This method was used, as opposed to a multi-layer tissue 
model, due to the dehydration of the tissue samples. Since it is difficult to know precisely the water content of each layer due 
to dehydration, a single-layer model allowed for easier manipulation of the contribution of water to the bulk absorption of the 
sample. The linear absorption coefficient of the sample defines the energy transfered from the optical source to the tissue. 

The BTEC model currently employs a single rate-process model of thermal injury. User-defined parameters for damage 
rates and activation energies as a function of temperature can be programmed for each tissue type or layer probed. A number 
of damage integral values or temperature shift searches are available in order to estimate damage thresholds or to compare to 
experimental data. 

For these simulations, the optical properties calculated based on equations for the absorption of a generic tissue were used 
in combination with the scattering coefficients measured by Salomatina et al.^^ Since absorption coefficient measurements 
for porcine skin for 1100-1550 nm are currently not found in literature, these coefficients were calculated using the equation 
for the net absorption coefficient of tissue proposed by Jacques.The net wavelength dependent absorption coefficient can 
be calculated from 

~ ^^f^a.oxy -f Jl(l S^^a.deoxy~^ ^ l^a.water ^ ^l^a.melanosome F 2.SC/,//,£/;/// (1) 

where S is the hemoglobin oxygen saturation of mixed arterio-venous vasculature, B is the average blood volume fraction, W 
is the water volume fraction, Bili is the bilirubin concentration, j3C is the j3-carotene concentration, F is the fat content, and 
M is the melanosome volume fraction. Due to their relatively low volume fraction and lack of literature values for absorption 
coefficients in the wavelength range used here, the contributions of j3-carotene, fat, hemoglobin, and bilirubin were not used 
in the calculation of the net absorption coefficient. 

Since the tissue samples used in this study were highly pigmented, thermal contributions due to melanin absorption must 
be taken into consideration. The wavelength-dependent absorption coefficient for melanin, Ba.meianosome^ was calculated using 

Ba.melanosome (519cn7 ) f 

where m is a power factor that is reported in literature to be approximately three.The bulk tissue sample was assumed to 
have a melanin volume fraction of 2.5% based upon previous reported values.Additionally, the differences in water content 
from sample to sample were accounted for using Eq 1 . Here, we allowed the water content to vary from 70% volume fraction 
to 35% volume fraction. Tissue perfusion was not incorporated into these simulations, since we had isolated non-perfused 
tissue. Damage thresholds were calculated for each wavelength by searching for a laser power that produced a prescribed 
value damage integral at a 100-/im depth in the tissue using the Arrhenius damage model.^*' 

Results 

The results of the experimentally-measured tissue heating are shown in Fig. 2. Temperature rises were calculated by subtract¬ 
ing the extracted maximum average temperature inside a 40-pixel (372 jj.m) diameter circular region of interest (ROI) centered 
at the middle of the laser spot from the average temperature in the same ROI prior to laser exposure. Fig. 1(B) shows a sample 
thermal image acquired for a 1425-nm laser exposure, taken shortly after the shutter opened. The black circle corresponds to 
the 40 pixel diameter ROI used to calculate the temperature rise. Extractions of temperature values from the ROI in the thermal 
images was done using a custom-written LabVIEW (National Instruments) virtual instrument (VI). A MATLAB script was 
used to compute the AT values for each exposure. 

The maximum AT was experimentally measured to be 31.8 ± 6.2 °C at 1425 nm. The significant deviation here can be 
attributed to changes in water volume fraction (VF/z^o) from sample to sample. The location of this peak corresponds to the 
strong water absorption band located in this wavelength region. Heating begins to fall off for wavelengths longer than 1450 
nm, as is expected based on water absorption and the continued decrease in melanin absorption. The minimum temperature 
rise measured was 1.9 ± 0.1 °C at 1100 nm. A slight increase in the thermal response is seen from 1125 nm to 1200 nm, 
which is expected based on the small increase in water absorption in this band. 

The results of the BTEC simulations are shown in Eig. 3a. Modeling results assumed a 36°C surface temperature for the 
tissue prior to laser exposure. The absorption coefficients used for the simulations were obtained by scaling the contribution 
from water absorption based on estimated changes in its volume fraction. 

The simulation results for VFh^o of 65% to 49% fit within the error bars measured in the experimental data. The tempera¬ 
ture rise for the 56% was found to be 34.5°C at 1425 nm. Simulations show a 37.6°C and 31.3°C temperature rise for 

65% VFhtO ^nd 49% VF/^^o respectively. The trends of the thermal response curves shown in Eig 3(A) match well with the 
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Figure 2. Experimentally-measured thermal response of the tissue samples to incident laser radiation. For all wavelengths, 
80 mW of average power was delivered to the sample. 
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Figure 3. Results of the BTEC simulations. (A) Simulated thermal response of tissue with varying (B) Estimated 

damage thresholds for select VFh^o- 
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experimental data, and show the importance of hydration when assessing damage limits, particularly when imaging biological 
samples and tissue with wavelengths in the IR. 

Fig. 3(B) shows the damage thresholds estimated for the absorption coefficient scaling that most closely matches the 
experimentally-measured temperature response. These trends follow what is expected based on the absorption coefficients 
of water and the measured thermal response. The use of 1100 nm through 1150 nm or 1250 nm through 1300 nm radiation 
allows for the most power to be delivered without causing photothermal tissue damage. These wavelengths allow for nearly a 
factor of 100 increase in laser power to be delivered to a tissue sample without causing damage as compared laser energies at 
1450 nm. 

Discussion 

As can be seen from the error bars in Fig 2, significant variation in the tissue heating was seen in the wavelength range of 1425 
to 1550 nm. Here, water is the predominate absorber, so slight changes in water content from sample to sample will result in 
the largest discrepancy. As tissue samples were excised and stored for varying times in the incubator prior to imaging, this too 
could result in significant changes in hydration. Slight variations in the relative humidity inside the acrylic box where the tissue 
samples were placed during imaging could also impact the measured thermal response. The data presented here highlight the 
importance of water content in a tissue sample with regard to the potential for photothermal damage. Measurement of the 
thermal response in an in vivo model would provide a more realistic assessment of damage thresholds with normal tissue 
hydration. Nonetheless, these data highlight the importance of hydration when assessing thermal damage related to laser 
exposure. 

Future modeling efforts should examine the effects of blood flow in addition to optical coefficients (/!„, /Xj, g) to explore 
how these parameters shape the trend in damage thresholds for photothermal damage. Additionally, further work is needed to 
provide accurate knowledge of how the optical properties of each tissue layer changes as a function of water volume fraction. 
With the increased prevalence of laser sources that produce output further into the infrared, such studies are critical to ensuring 
safe operation across all realms of use. 
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